UNCLASSIFIED/UNLIMITED 


ORGANIZATION 


Waveform  Diversity:  Past,  Present,  and  Future 

A.  De  Maio  (SIEEE),  A.  Farina  (FREng,  FIET,  FffiEE) 


Abstract 

Waveform  diversity  indicates  the  ability  to  adapt  and  diversify  dynamically  the  waveform  to  the 
operating  environment  in  order  to  achieve  a  performance  gain  over  non-adaptive  systems.  This  technique 
can  allow  one  or  more  sensors  to  automatically  change  operating  parameters  such  as  frequency,  pulse 
repetition  time,  transmit  pattern,  modulation,  etc.  The  present  lecture  starts  with  an  overview  concerning 
the  role  of  the  waveform  diversity  in  history,  mathematics,  and  music  from  the  epoch  of  Pythagoras, 
continuing  with  the  studies  of  Galileo,  Fourier,  and  Maxwell.  Examples  of  waveform  diversity  in  nature, 
such  as  the  bath  sonar  signal,  the  sounds  of  whales,  and  the  cosmic  microwave  background  radiation 
are  presented1.  A  tutorial  introduction  to  the  concept  of  ambiguity  function,  its  relevant  properties,  and 
its  role  as  an  instrument  to  quantify  the  quality  of  a  waveform,  follows.  Precisely,  after  a  short  review  of 
the  most  common  radar  signals  and  their  ambiguity  functions,  the  effects  of  a  possible  signal  coding  is 
thoroughly  described.  Amplitude,  phase,  and  frequency  codes  are  considered,  even  if  a  special  attention 
is  deserved  to  the  class  of  frequency  coded  waveforms  through  a  Costas  sequence. 

Keywords.  Ambiguity  Function,  Radar  Coding,  Coherent  Train  of  Diverse  Pulses. 

I.  Introduction 

The  waveform  exploited  by  the  radar  is  responsible  of  resolution,  accuracy,  and  ambiguity  of 
the  target  range  and  radial  velocity  measurements.  While  range  is  associated  with  the  delay  of 
the  received  signal,  radial  velocity  depends  on  the  Doppler  frequency  shift. 

Waveform  design  algorithms  usually  anticipated  their  implementation  by  many  years,  due 

to  complexity  and  hardware  limitations  [1].  For  instance,  the  concept  of  pulse  compression, 
developed  during  the  World  War  II,  gained  renewed  interest  only  when  high-power  Klystrons 
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became  available  [2],  In  other  words,  what  seems  unpractical  today,  may  not  be  definitely  ruled 
out  in  the  near  future.  The  lack  of  signal  coherence,  which  precluded  the  application  of  signal 
compression  during  the  World  War  II  is  today  easy.  Maybe,  the  linear  power  amplifiers,  required 
to  implement  amplitude  modulated  radar  signals,  will  not  represent  a  technological  limitation 
tomorrow. 

If  a  matched  filter  is  used  at  the  receiver,  the  ambiguity  function  represents  a  suitable  tool 
to  study  the  response  of  the  filter  in  two  dimensions:  delay  and  Doppler.  The  constant  volume 
underneath  the  squared  ambiguity  function  involves  some  trade-offs  in  signal  design.  Precisely, 
a  narrow  response  in  one  dimension  is  accompanied  by  a  poor  response  in  the  other  dimension 
or  by  additional  ambiguous  peaks.  Moreover,  if  we  prefer  ambiguous  peaks  to  be  well  spaced 
in  delay,  we  have  to  accept  them  closely  spaced  in  Doppler  (and  vice  versa).  If  we  want  a  good 
Doppler  resolution,  we  need  long  coherent  signal  durations. 

Several  signals  are  used  for  different  radar  applications  and  systems.  Modern  pulsed  radars 
generally  use  pulse  compression  waveforms  characterized  by  high  pulse  energy  (with  no  increase 
in  peak  power)  and  large  pulse  bandwidth.  As  a  consequence,  they  provide  high  range  resolution 
without  sacrificing  maximum  range  which  depends  on  the  pulse  energy. 

Unfortunately,  there  are  not  easily-handled  mathematical  techniques  to  calculate  a  signal  with 
a  prescribed  ambiguity  function.  It  follows  that  the  design  of  a  radar  signal  with  desirable 
characteristics  of  the  ambiguity  function  is  mainly  based  on  the  designer’s  prior  knowledge  of 
radar  signatures  as  well  as  on  “trial  and  check  procedures”. 

In  this  lecture,  we  first  present  (Section  II)  the  mathematical  definition  of  the  ambiguity 
function  and  describe  its  relevant  properties.  Then,  we  will  explore,  in  Section  III,  the  ambiguity 
function  of  some  basic  radar  signals:  single-frequency  rectangular  pulse,  Linear  Frequency 
Modulated  (LFM)  pulse,  and  coherent  pulse  train.  Hence,  in  Section  IV,  the  conflicts  in  designing 
suitable  waveforms  for  different  applications  are  discussed:  radar  coding  is  presented  as  a  suitable 
mean  to  achieve  ambiguity  function  shaping.  Several  techniques  based  on  frequency  and  phase 
coding  are  presented;  the  ultimate  goal  is  to  segregate  the  volume  of  the  ambiguity  function 
in  regions  of  the  delay-Doppler  plane  where  it  ceases  to  be  a  practical  embarrassment  [3],  In 
Section  V,  the  merits  and  the  drawbacks  concerning  the  use  of  coherent  trains  of  diverse  pulses 
is  addressed.  Finally  concluding  remarks  are  drawn  in  Section  VI. 
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II.  Ambiguity  Function:  Definition  and  Properties 


This  function  was  introduced  in  signal  analysis  by  Ville  [4]  and  in  the  radar  context  by 
Woodward  [3].  However,  it  was  known  in  thermodynamic,  since  1932,  due  to  Eugene  Wigner 
(Nobel  prize)  who  studied  quantum  corrections  to  classical  statistical  mechanics  [5], 

The  ambiguity  function  of  a  signal  whose  complex  envelope  is  denoted  by  u(t )  is  defined  as 


\X(t,v)  | 


u(t)u*(t  +  r)  exp(j2nut)dt 


(1) 


where  (•)*  represents  the  conjugate  operator,  j  =  \/—l,  |  - 1  is  the  modulus  of  a  complex  number, 
r  and  v  are  the  incremental  delay  and  Doppler  frequency  shift  respectively.  Otherwise  stated, 
it  is  the  modulus  of  a  matched  filter  output  when  the  input  is  a  Doppler  shifted  version  of  the 
original  signal  to  which  the  filter  is  actually  matched.  It  follows  that  |X(0,  0)  |  coincides  with  the 
output  when  the  input  signal  is  matched  to  the  nominal  delay  and  Doppler  of  the  filter;  non-zero 
values  of  r  and  v  indicate  a  target  from  other  range  and/or  velocity. 

Assuming  that  u(t)  has  unitary  energy,  |A(r,  v)\  complies  with  the  following  four  relevant 
properties. 

1)  Maximum  Value  Property. 


\X(t,u)\<\X(0,  0)|  =  1, 


(2) 


the  maximum  value  of  the  ambiguity  function  is  reached  for  (r,  u) 
to  1. 

2)  Unitary  Volume  Property. 


| X(t,  v)\2  drdu  =  1 , 


'  —  OO  J  —  OO 


the  volume  underneath  the  squared  ambiguity  function  is  unitary. 

3)  Symmetry. 


\X(t,p)\  =  \X(-t,-p)\  , 


(0, 0)  and  is  equal 


(3) 


(4) 


the  ambiguity  function  shares  a  symmetry  property  about  the  origin. 

4)  Linear  Frequency  Modulation  (LFM)  Property. 

If  \X(t,  v)\  is  the  ambiguity  function  corresponding  to  u(t),  then  \X(r,u  —  kr)\  is  the 
ambiguity  function  of  u(t)  exp(jirkt2). 
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A  more  concise  way  of  representing  the  ambiguity  function  consists  of  examining  the  one¬ 
dimensional  zero-delay  and  zero-Doppler  cuts.  The  cut  of  | X(t,  v)  along  the  delay  axis  is 


|X(r,0)|  = 


u(t)u*(t  +  r)dt 


where  R(r)  is  the  autocorrelation  function  of  u(t).  The  cut  along  the  Doppler  axis  is 


l*(o,OI  = 


\u(t)\2  exp(j2TTut)dt 


(5) 

(6) 


which  is  independent  of  any  phase  or  frequency  modulation  of  the  input  signal.  Further  interesting 
properties  of  the  ambiguity  function  can  be  found  in  [6],  Finally,  in  [1],  the  concept  of  periodic 
ambiguity  function  is  presented  and  its  connection  with  (1)  is  discussed. 


III.  Ambiguity  Function  of  Basic  Radar  Signals 

In  this  section,  we  present  the  ambiguity  function  of  some  basic  signals  (single  frequency 
rectangular  pulse,  LFM  pulse,  and  coherent  pulse  train)  [7,  eh.  8]  and  discuss  their  suitability 
for  radar  applications. 


A.  Rectangular  Pulse 


The  rectangular  pulse  of  length  tp  and  unitary  energy  is  given  by1 


.  ,  1  ft 

u(t)  =  ^=rect  — 


and  the  corresponding  ambiguity  function  is 


I  X(r,u)\  =  { 


1  _  \ll\  sinAAC1  “  MApM 


tp  J  ntp{  1  -  \r\/tp)u 
0  elsewhere 


T 


1 - )  sine  [tp(  1  — 


MApH 


I T I  t 

J  1 1  I  _  Lp  5 


(7) 


In  Figures  la-lc,  (7)  is  plotted  together  with  the  contours  and  the  cuts  along  the  delay  and 
Doppler  axes.  Notice  that  (7)  is  limited  to  an  infinite  strip  whose  size  on  the  delay  axis  is  2 tp. 
As  to  the  cut  at  r  =  0,  it  exhibits  the  first  nulls  at  vnun  =  ±7-  and,  since  the  sinc(-)  function 
has  a  peak  sidelobe  at  —13.5  dB,  the  practical  extension  of  the  ambiguity  function  along  the 
Doppler  axis  can  be  considered  2 /tp. 
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Figure  la:  Ambiguity  function  of  a  constant  frequency  rectangular  pulse  of  length  tp. 


In  general,  the  square  pulse  is  not  a  desirable  waveform  from  a  pulse  compression  standpoint, 
because  the  autocorrelation  function  is  too  wide  in  time,  making  it  difficult  to  discern  multiple 


lrThe  function  rect(x)  is  equal  to  1,  if  |x|  <  1/2,  and  is  equal  to  0  elsewhere.  The  function  sinc(:c)  is  defined  as  sinc(x)  = 
sin(7rx) 

7 TX 
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Figure  lc:  Ambiguity  function  of  a  constant  frequency  rectangular  pulse  of  length  tp.  a)  Zero-Doppler  cut.  b) 

Zero-delay  cut. 


overlapping  targets. 


B.  LFM  Pulse 


The  LFM  pulse  or  chirp  is  commonly  used  in  radar  and  sonar  applications.  It  has  the  advantage 
of  greater  bandwidth  while  keeping  the  pulse  duration  short  and  the  envelope  constant.  The 
complex  envelope  of  a  LFM  pulse,  with  instantaneous  frequency  f(t )  =  kt,  is 


u(t)  = 


rect  (  —  )  exp(jirkt2) , 


t. 


and  the  corresponding  ambiguity  function  is  given  by 


\X{t,v)\  = 


1  - 


T 


sine  [tp(l  —  \r\/tp){v  —  kr)\ 


\r\  <  tp , 


(8) 


0  elsewhere 

\ 

In  Figures  2a-2c,  (8)  is  plotted  together  with  the  contours  and  the  cuts  along  the  delay  and 
Doppler  axes. 

Notice  that  the  cut  along  the  Doppler  axis  (r  =  0)  is  the  same  as  in  Figure  lc-b.  On  the 
contrary,  the  cut  along  the  Delay  axis  (y  =  0)  is  deeply  different  from  Figure  lc-a:  if  kt2p  = 
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Figure  2a:  Ambiguity  function  of  a  LFM  rectangular  pulse  of  length  tp  and  with  kt*  =  10. 


tPAf  4  (A /  is  the  total  frequency  deviation),  it  exhibits  the  first  nulls  at 

1  1 

Tnull  i  7  ,  i  A  7  • 
ktp  A  / 

This  means  that  the  range  window  has  been  compressed  by  a  factor  D  =  tpAf,  which  is  usually 
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Figure  2c:  Ambiguity  function  of  a  LFM  rectangular  pulse  of  length  tp  and  with  ktp  =  10.  a)  Zero-Doppler  cut. 

b)  Zero-delay  cut. 


referred  to  as  compression  ratio.  Notice  also  that  the  ambiguity  function  volume  is  mainly 
concentrated  on  a  diagonal  ridge. 

Slight  Doppler  mismatches  for  the  LFM  pulse  do  not  change  the  general  shape  of  the  pulse 
and  reduce  the  amplitude  very  little,  but  they  appear  to  shift  the  pulse  in  time.  Thus,  an 
uncompensated  Doppler  shift  changes  the  target’s  apparent  range;  this  phenomenon  is  called 
range-Doppler  coupling. 

Finally,  we  just  mention  that  non-linear  FM  pulses  can  be  conceived  (see  for  instance  [8]  and 

[1]). 

C.  Coherent  Pulse  Train 

The  complex  envelope  of  a  coherent  pulse  train,  composed  by  N  equally  spaced  pulses,  can 
be  written  as 

1  N 

u(t)  =  -j=^2un(t-  (n-  1)Tr)  (9) 

v  n= 1 

where  TR  is  the  pulse  repetition  period  and  un(t)  is  the  complex  envelope  of  the  n-th  unitary 
energy  pulse.  Assuming  that  the  pulse  train  is  uniform  (i.e.  un(t)  =  uc(t),  n  =  1 , ,N)  and 
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that  the  separation  between  pulses  TR/2  is  greater  than  the  pulse  duration  tp,  the  ambiguity 
function  of  (9)  can  be  expressed  as 

i  N~ 1 

I X(T’^l  =  iV  ^2  \Xc(r  -pTR,u)\ 

P=-(N- 1) 

where  | Xc(r,  v)\  is  the  ambiguity  function  of  ucif). 

In  Figure  3,  we  assume  single-frequency  rectangular  pulses,  N  =  6,  TR  =  5 tp  and  plot  (10) 
in  the  range-Doppler  domain2.  Due  to  its  shape  (10)  is  often  referred  to  as  bed  of  nails.  The 
zero-Doppler  cut  shows  that  there  are  multiple  triangular  windows:  the  separation  between  two 
consecutive  peaks  is  equal  to  the  pulse  repetition  period  TR.  Moreover,  all  the  triangular  windows 
have  the  same  width  2 tp,  but  their  height  decreases  as  the  distance  from  the  origin  increases. 

As  to  the  cut  for  r  =  0,  there  are  multiple  peaks  spaced  apart  1  /TR  and  N  —  2  smaller 
sidelobes  between  them.  The  first  nulls  occur  at  v  =  ±1  /NTR,  namely  the  width  of  the  main 
peak  (in  Doppler)  is  ruled  by  the  length  of  the  Coherent  Processing  Interval  (CPI). 


sin[7ri/(iV  -  \p\)TR] 


sm(nuTR) 


(10) 


Figure  3:  Ambiguity  function  of  a  coherent  train  of  uniform  pulses  with  N  =  6,  pulse  length  tp ,  and  pulse 

repetition  period  Tr  =  5 tp. 


2In  the  following,  the  MATLAB  toolbox  of  [9]  is  used  to  plot  the  ambiguity  functions. 
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IV.  Coded  Radar  Signals 

The  ambiguity  function  of  a  coherent  pulse  train  allows  a  main  peak  narrow  both  in  range  and 
in  Doppler,  but  exhibits  some  peaks  with  almost  the  same  amplitude  as  the  main  peak.  These 
might  be  deleterious  and  can  lead  to  range/Doppler  ambiguities  very  difficult  to  resolve. 

If  we  wish  to  maintain  a  very  narrow  main  peak  but  cannot  accept  the  additional  peaks  typical 
of  the  bed  of  nails,  we  can  spread  the  volume  in  a  low  but  wide  pedestal  around  the  main  peak. 
This  kind  of  ambiguity  function  is  referred  to  as  thumbtack  shape  and  can  be  obtained  considering 
coded  radar  signals. 

A.  Frequency  Coding:  Costas  Sequences 

The  complex  envelope  of  a  frequency  coded  pulse  of  length  tp  can  be  written  as 

1  N 

u{t)  =  ~7Wr^2Un^~  1)^)>  (u) 

v  - Ntb  n=1 

where 

f  exp(j27r  fnt)  0  <t  <tb 

Unit)  =  \  (12) 

I  0  elsewhere 

4  is  the  length  of  each  subpulse  (time-slot  duration,  Nf  =  tp),  the  frequency  shift  in  the  n-th 
time  slot  is  fn  =  an/tb,  while  the  hopping  (coding)  sequence  is 

®1>  •  •  •  >  diV  )  un  £  {0, . . . ,  N  1)  • 

The  frequency  history  of  (11)  can  be  represented  through  the  coding  matrix  (Table  I)  where  the 
horizontal  axis,  representing  time,  is  divided  in  N  time-slots  of  length  4  and  the  vertical  axis 
is  used  to  represent  equally  spaced  frequencies.  The  {h,  k)- th  entry  of  the  binary  matrix  can 
assume  only  two  values:  1  if  the  h- th  frequency  is  transmitted  in  the  A  -th  time  slot,  0  elsewhere. 
Obviously,  there  is  only  a  1  per  column. 
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0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

Frequency 


Time 


Table  I:  Binary  matrix  representation  of  frequency  coding. 


The  corresponding  ambiguity  function  can  be  evaluated  through  the  expression 


\X(t,v)\  = 


1 

N 


N 


"IP  exp(j2n(m  -  1  )utb) 


m=  1 


N 


$mm(T,  U)  +  ^  $mn(T  -  (m  -  n)tb,  u) 

n=l,m/n 


(13) 


where 


\T\ 

tb 


sine (amn)  exp (-j(3mn  -  j2irfnT),  |r|  <  tb  , 


(14) 


0  elsewhere 


and 


®-mn  (fm  fn  V)  (tb  |^"|) 
fimn  7T  (fm  fn  ^)  (fft  T  7”)  . 

Slightly  different  codes  can  strongly  affect  the  ambiguity  function  of  the  signal;  hence  it  is  of 
interest  to  present  a  methodology  which  roughly  predicts  the  ambiguity  shape.  Such  a  technique 
is  based  on  the  observation  that  the  cross  correlation  between  signals  at  different  frequencies 
approaches  zero  when  the  frequency  difference  is  large  with  respect  to  the  inverse  of  the  signal 
duration  (or  equal  to  multiples  of  that  inverse).  The  prediction  is  possible  overlaying  a  copy  of 
the  binary  matrix  on  itself,  and  then  shifting  one  relative  to  the  other  according  to  the  desired 
delay  (horizontal  shifts)  and  Doppler  (vertical  shifts).  A  coincidence  of  two  elements  in  the 
matrix  denotes  a  peak  of  amplitude  one  in  the  predicted  ambiguity  function,  two  coincidences 
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Figure  4a:  Ambiguity  function  of  a  Costas  signal  with  coding  sequence  (3,  6,  0,  5,  4,  1,  2). 


a  peak  of  amplitude  2,  and  so  on.  The  maximum  number  of  coincidences  is  the  number  of 
frequencies  (N  in  our  example),  and  can  be  reached  only  in  the  origin  (zero  delay  and  zero 
Doppler).  Normalizing  the  maximum  peak  at  1,  we  can  assume  a  coincidence  equal  to  a  peak 
of  amplitude 

Definition  1.  A  coding  sequence  is  a  Costas  code  [10]  if  all  the  non-zero  shifts  of  the  binary 
matrix  do  not  lead  to  more  than  one  coincidence. 

In  Figure  4a,  we  plot  the  ambiguity  function  of  a  Costas  coded  pulse  with  N  —  7  and  coding 
sequence  (oi, . . .  ,ayv)  =  (3,  6,  0,  5,  4,  1,  2).  The  thumbtack  nature  of  the  ambiguity  function 
is  clearly  evident.  Moreover,  in  Figure  4b,  we  plot  the  contours  of  (13)  for  | X(t,  p)  =  0.125. 
Notice  that  there  is  a  similarity  between  Figure  4b  and  the  sidelobe  matrix  [1,  p.  77]  of  the 
coding  sequence  (Table  II).  In  Figure  4c,  we  plot  the  autocorrelation  function  (zero-Doppler 
cut):  as  expected,  there  are  nulls  at  r  =  kt &.  Moreover,  according  to  Property  4  of  the  ambiguity 
function,  we  do  not  need  to  plot  the  zero-delay  cut,  since  it  does  not  depends  on  the  frequency 
modulation,  but  only  on  the  magnitude  of  the  unmodulated  pulse. 
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Figure  4b:  Ambiguity  function  contours  at  0.125  of  a  Costas  signal  with  coding  sequence  (3,  6,  0,  5,  4,  1,  2). 
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—6 

-5 

-4 

-3 

-2 

-1 
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+1 

+2 

+3 

+4 

+5 

+6 

Delay 


Table  II:  Sidelobe  matrix  of  a  Costas  signal  with  coding  sequence  (3,  6,  0,  5,  4,  1,  2). 


Unfortunately,  it  does  not  exist  a  constructive  procedure  to  determine  all  the  possible  Costas 
sequences  of  a  fixed  length,  nor  how  many  they  are.  To  circumvent  this  drawback,  two  approaches 
can  be  followed. 

•  An  exhaustive  search  among  all  the  N\  possible  sequences  of  length  N; 

•  A  constructive  procedure  to  determine  a  subclass  of  particular  Costas  sequences. 

The  first  method  needs  grid  computing,  and  cannot  provide  very  long  sequences.  For  example, 
using  a  grid  of  more  than  700  processors,  a  complete  search  of  Costas  array  of  length  30 
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Figure  4c:  Autocorrelation  function  of  a  Costas  signal  with  coding  sequence  (3,  6,  0,  5,  4,  1,  2). 


requires  more  than  4  years  [11].  Actually,  the  public  database  of  Costas  arrays  contains  all  the 
sequences  starting  from  the  unique  sequence  of  length  1,  up  to  the  204  sequences  of  length  27 
[1 1],  Moreover,  the  ratio  between  the  number  of  length  N  Costas  sequences  and  the  N\  possible 
sequences,  decreases  very  quickly  [11],  namely  it  is  even  more  difficult  to  find  a  Costas  sequence 
increasing  the  length. 

A  different  construction  technique  is  based  on  the  theory  of  Galois  finite  fields3.  Starting 
from  a  primitive  element  of  GF(N),  i.e.  an  element  of  the  field  that  can  generate  all  the  others 
elements  but  for  0,  it  is  possible  to  conceive  several  procedures  to  construct  a  Costas  sequence. 
The  most  used  techniques  are  the  Welch  1,  the  Welch  2,  the  Golomb  2,  the  Lempel  2,  and  the 
Taylor  4.  Let  us  now  illustrate  how  the  Welch  1  procedure  can  be  implemented. 

Choose  a  length  N  >  3  such  that  N  =  p  —  1,  where  p  is  a  prime  number.  Find  a  primitive 
element  a  of  GF(p).  Numbering  the  columns  of  the  array  in  Table  III  with  k  =  0, 1,2, . .  .,p- 2 
and  the  rows  with  h  =  1, 2, . . .  ,p  —  1,  we  put  a  1  in  position  (h,  k )  if  and  only  if  h  =  ak  (mod 
p).  For  example,  let  us  consider  N  =  4,  so  p  =  5.  A  primitive  element  of  GF( 5)  is  2,  since 
the  elements  {1, 2, 3, 4}  can  be  obtained  as  {2°,  21, 23, 22}  (mod  p).  Now,  we  can  construct  the 

3In  the  following,  a  Galois  field  containing  the  elements  from  0  to  N  —  1  will  be  denote  by  GF(N). 
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Frequency 


0 

1 

2 

CO 

1 

1 

0 

0 

0 

2 

0 

1 

0 

0 

CO 

0 

0 

0 

1 

4 

0 

0 

1 

0 

Time 


Table  III:  Welch  I  construction  matrix. 


which  provides  the  Costas  code  (1,  2,  4,  3). 

Finally,  in  Table  IV,  we  present  a  short  list  of  procedures,  which  can  be  used  to  obtain  Costas 
sequences  whose  length  ranges  between  3  and  30  [12]. 


Order 

Working  constructions 

Order 

Working  constructions 

— 

— 

16 

Ti,  Wi,  W3,  L3,  G3 

— 

— 

17 

To,  W2,  L2,  G2 

3 

T\,  W2,  L2,  G2 

18 

Wi 

4 

Ti,  Wi,  G3,  G4,  G5 

19 

Wo 

5 

T0,W2iL2,T4,G2,G3,G4 

20 

g3 

6 

t4,w3,l2,g2,g3 

21 

w2,l2,g2 

7 

w0,l2,t4,g2 

22 

tuwug3 

8 

t1,w3,l3ig3 

23 

To,  L2,  G2 

9 

w2,l2,g2 

24 

g3 

10 

t2,w1,w3,l3,g3 

25 

l2,  g2 

11 

t0,w2,l2,g2 

26 

w3,l2,g2 

12 

w4,g4 

27 

w2,l2,t4,g2 

13 

w0,g3 

28 

Ti,  W\,  G3,  G4 

14 

l2,g2,  g3 

29 

To,W2,L2,G2,G3 

15 

w2,l2,t4,g2 

30 

w3,l2,g2 

Table  IV 


Constructions  that  successfully  produce  Costas  arrays  of  order  <  30. 
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Figure  5:  Ambiguity  function  of  a  Costas  sequence  of  length  40. 


The  ambiguity  function  behavior  of  the  Costas  coding  is  more  pronounced  for  long  sequences. 
In  Figure  5,  we  plot  the  ambiguity  function  of  a  Costas  sequence  of  length  40  [1].  This  last 
example  of  a  relatively  long  Costas  signal  indicates  that  the  longer  the  sequence,  the  closer  its 
ambiguity  function  to  the  thumbtack  shape.  Regarding  to  the  zero-Doppler  cut,  the  sidelobes 
are  usually  below  —26  dB.  However,  the  near  sidelobes  (h/N  <  r  <  t\,)  are  higher,  decaying 
from  —13.7  dB  in  a  manner  typical  of  the  autocorrelation  function  sidelobes  of  a  signal  with 
a  rectangular  spectrum.  Indeed,  the  spectrum  of  our  relatively  long  Costas  signal  is  nearly 
rectangular. 

B.  Frequency  Coding:  Pushing  Sequences 

This  class  of  frequency  coded  signals  has  been  introduced  by  Chang  and  Bell  in  [13], 

Definition  2.  For  the  ambiguity  function  of  a  frequency  coded  waveform,  a  clear  area  of  size 
s  is  a  connected  area  centered  at  the  origin  of  the  (r,  i^-plane,  where  |r|  <  sti,  and  \v\  <  s/tb  , 
such  that  no  sidelobe  peaks  (of  height  greater  than  1  /N)  are  present  in  this  area. 

A  frequency-coding  sequence  having  the  ambiguity  function  with  a  clear  area  of  size  s  is 
called  a  pushing  sequence  with  power  s,  where  s  >  1. 

Of  particular  relevance  are  the  codes  which  are  both  Costas  and  pushing  (Costas-pushing 
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Figure  6a:  Ambiguity  function  of  a  Costas-pushing  sequence  of  length  5. 


sequences).  They  share  some  interesting  symmetry  properties  (Group  D{  Dihedral  Symmetry 
Property)  [13]  which  permit  to  derive  more  pushing-Costas  sequences  from  a  given  one.  More¬ 
over  a  Costas  sequence,  derived  from  the  Taylor  4  (T±)  construction,  is  also  a  pushing  sequence 
of  power  1. 

In  Figures  6a-6b  we  plot  the  ambiguity  function  and  the  contours  at  \X(t,v)\  =  0.1667  of 
a  Costas-pushing  coded  pulse  with  N  =  5  and  coding  sequence  (ai, . . . ,  a/v)  =  (1,  3,  0,  4,  2). 
The  clear  area  and  the  thumbtack  nature  of  the  ambiguity  function  can  be  directly  seen  from 
the  figures. 

C.  Phase-Coded  Radar  Signals 

A  pulse  of  duration  tp  is  divided  into  N  bits  (chips)  of  identical  duration  f6  =  tp/N  and  each 
of  them  is  coded  with  a  different  pulse  value.  The  complex  envelope  of  the  phase-coded  pulse 
is  given  by 


(15) 


where  un  =  exp (j(f>n)  and  the  set  of  N  phases  02, . . . ,  4>n}  is  the  phase  code  associated  with 
u(t).  Criteria  for  selecting  a  specific  code  are  the  resolution  properties  of  the  resulting  waveform 
(shape  or  ambiguity  function),  frequency  spectrum  and  the  ease  through  which  the  system  can  be 
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Figure  6b:  Ambiguity  function  contours  at  0.1667  of  a  Costas-pushing  signal  with  code  sequence  (1,  3,  0,  4,  2). 


implemented.  However,  finding  a  code  which  leads  to  a  predetermined  range-Doppler  resolution 
is  very  complicated;  choosing  a  code  ensuring  a  good  autocorrelation  function  is  easier. 

The  autocorrelation  function  of  a  phase-coded  pulse  is  a  continuous  function  of  the  delay  r. 
The  properties  of  the  autocorrelation  function  should  be  examined,  in  general,  for  all  —tp< 
T  <  tp.  It  can  be  shown  [1]  that  it  is  sufficient  to  calculate  the  correlation  function  at  integer 
multiples  of  the  bit  duration.  The  other  can  obtained  connecting  the  values  at  r  =  kR  using 
straight  lines  in  the  complex  plane.  Thus,  the  optimization  of  the  continuous  autocorrelation 
peaks  is  simplified  to  the  minimization  of  the  discrete  correlation  function  \Rk\  values  [1]. 

1 )  Barker  codes:  The  Barker  codes  [14]  were  designed  as  the  sets  of  N  binary  pulses  yielding 
a  peak-to-peak  sidelobe  ratio  equal  to  N .  In  [14]  and  [15]  all  the  known  binary  sequences 
complying  with  this  property  are  reported.  It  has  been  shown  that  no  Barker  codes  exist  for 
13  <  N  <  1.898.884  and  for  all  odd  N  >  13  [16].  In  Figure  7a  the  autocorrelation  function  for 
a  Barker  code  of  length  iV  =  13  is  shown.  In  Table  V,  all  the  known  binary  Barker  codes  are 
reported. 

One  of  the  main  drawbacks  of  the  Barker  codes  is  that  they  are  optimized  only  for  the 
autocorrelation  function,  i.e.  the  chosen  peak-to-peak  sidelobe  ratio  is  valid  only  for  the  zero- 
Doppler  cut.  Thus,  if  the  target  return  is  Doppler  shifted  the  sidelobe  peaks  are  greater  than  the 
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Figure  7a:  Autocorrelation  function  of  a  phase-coded  pulse  using  the  13-element  Barker  code. 


Code  length 

Code 

2 

11  or  10 

3 

110 

4 

1110  or  1101 

5 

11101 

7 

1110010 

11 

11100010010 

13 

1111100110101 

Table  V 


All  known  binary  Barker  codes. 


expected.  Otherwise  stated,  Barker  codes  exhibit  a  low  Doppler  tolerance.  Figure  7b  shows  the 
ambiguity  function  of  a  13-element  Barker  code;  notice  that  the  sidelobe  peaks  are  higher  than 
1/N  for  many  values  of  the  Doppler  shift  and  the  main  peak  is  lower  than  1  for  those  Doppler 
values. 
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Figure  7b:  Ambiguity  function  for  a  13-element  Barker  code. 


2)  Polyphase  Barker  Codes:  The  Barker  codes  are  limited  by  the  binary  assumption.  Allowing 
any  phase  value  can  lead  to  lower  sidelobes.  In  this  way,  a  polyphase  code  is  realized.  The 
polyphase  iV-sequence  with  minimal  peak-to-sidelobe  ratio  excluding  the  outermost  sidelobe 
(which  is  always  1  /N,  both  for  binary  and  polyphase  codes)  is  called  generalized  Barker 
sequence  or  polyphase  Barker  sequence.  Actually,  systematic  methods  to  produce  polyphase 
Barker  sequences  are  not  yet  found;  using  numerical  optimization  techniques  allow  to  search 
such  codes  without  restrictions  on  the  values  of  the  sequence  phases.  Figure  8  shows  the 
autocorrelation  function  magnitude  for  the  13-element  polyphase  Barker  code.  Notice  that  the 
peak  sidelode,  whose  magnitude  is  1/N,  is  located  at  r  =  ±12/4. 

3)  Frank  Codes  and  extensions:  The  Frank  code  [17]  is  derived  from  the  phase  history  of  a 
linearly  frequency  stepped  pulse;  the  Frank  code  is  designed  to  ensure  low  sidelobe  peaks  even 
for  non  zero-Doppler  values.  However,  it  applies  only  for  a  square  code  length  (N  =  L2).  The 
elements  un  (1  <  n  <  N)  of  a  N  =  L2  Frank  code  are 

U(n-i )L+k  =  exp (j(j)n,k)  1  <  n  <  L  1  <k  <  L  (16) 
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Figure  8:  Autocorrelation  function  of  13-element  polyphase  Barker  code. 


where  ([)n±  =  2n(n  —  1  )(k  —  1  )/L.  The  phase  values  in  (16)  can  be  obtained  from  the  elements 
of  the  L  x  L  discrete  matrix 

"  0  0  0  ...  o’ 

0  1  2  ...  L- 1 

0  2  4  ...  2(L-1)  , 

_  0  (L  -  1)  2(L  —  1)  ...  (L-  l)2  _ 

concatenating  the  rows  and  multiplying  for  2i r/L;  finally,  the  phase  values  are  taken  mod  2n. 
The  Frank  code  has  two  important  properties  [1]: 

1)  it  is  perfect4; 

2)  the  autocorrelation  function  exhibits  relatively  low  sidelobes. 

Figure  9  shows  the  ambiguity  function  for  a  16-elements  Frank  code. 

4A  phase  code  having  zero  periodic  autocorrelation  sidelobes  is  called  perfect. 
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Figure  9:  Ambiguity  function  of  16-element  Frank  code. 


The  Frank  code  can  be  modified  into  the  PI,  P2  and  PX  codes,  with  the  dc  frequency5  term 
in  the  middle  of  the  pulse  instead  of  the  beginning;  also  the  modified  versions  can  be  applied 
only  for  square  code  length,  i.e.  N  =  L2.  The  PX  code  [18]  yields  the  same  aperiodic  peak 
sidelobe  as  the  Frank  code,  but  ensures  a  low  integrated  sidelobe  level.  The  elements  of  the  PX 
code  are  defined  as 


•S(n—  l)L+k  6Xp  (^j4>n,k) 


(17) 


where 


L  even 


L  odd . 


(18) 


In  the  same  way  as  for  the  Frank  code,  it  is  possible  to  derive  a  matrix  to  form  the  phase  values 
for  the  PX  code. 

P2  code  [19] -[20]  can  be  constructed  only  for  L  even  and  is  defined  exactly  as  the  PX  code. 
This  code  is  palindromic  as  it  exhibits  some  specific  symmetry  properties  [1].  PI  code  [19] -[20] 


5 


The  dc  frequency  term  is  in  correspondence  to  the  zero-phase  value  in  the  code. 
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elements  are  defined  using  (17)  where 

4>n,k  =  jr  [(n  -  l)L  +(k-  1)]  •  (19) 

The  PI  code,  unlike  PX  and  P2  codes,  is  perfect,  as  the  Frank  code.  The  ambiguity  function  of 
the  PI  code  for  odd  L  is  identical  to  that  of  the  Frank  code.  For  even  L,  the  ambiguity  functions 
of  P2  and  PX  codes  are  very  similar  to  the  one  of  the  PI  code  and  also  to  the  one  of  the  Frank 
code. 


D.  A  Table  with  Some  Common  Radar  Codes 

There  exist  further  coding  strategies,  different  from  those  presented  in  the  previous  sub¬ 
sections.  The  interested  reader  can  consult  the  excellent  book  of  Levanon  and  Mozeson  [1]  for 
a  more  complete  description.  Here,  we  just  provide  a  table  (Table  VI,  of  course  not  exhaustive) 
to  summarize  some  common  radar  codes  [21]-[35]. 
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Name 

Inventor 

Year 

Type 

Barker  Code 

R.H.  Barker 

1953 

Phase  Codes 

Complementary  Code 

M.J.E.  Golay 

1961 

Complementary  Phase  Codes 

Huffmann  Code 

D.A.Huffmann 

1962 

Phase  Code 

Frank  Code 

R.L.  Frank,  S.A.  Zadoff 

1962 

Chirplike  Phase  Codes 

Zadoff-Chu  Code 

S.A.  Zadoff 

1963 

Chirplike  Phase  Codes 

Gold  Code 

R.  Gold 

1967 

Phased  Code  based  on  binary  sequence 

Minimum  peaksidelobe  code 

J.  Lindner,  N.  Cohen  et  al. 

1975 

Phased  Code  based  on  binary  sequence 

Welti  Code 

R.  Sivaswamy 

1978 

Subcomplementary  Code 

PI  and  P2  Codes 

B.L.  Lewis,  F.F.  Kretschmer 

1981 

Phase  Codes 

Frank  Polyphase  Codes 

B.L.  Lewis,  F.F.  Kretschmer 

1983 

Polyphase  Codes 

Costas  Array 

J.P.  Costas 

1984 

Frequency  Codes 

Quadratic  congruential  coding 

J.R.  Bellegarda,  E.L.  Titlebaum 

1988 

Frequency  Codes 

Polyphase  Barker  Codes 

L.  Bomer,  M.  Autweiler 

1989 

Phase  Codes 

Generalized  P4  Code 

F.F.  Kretshmer,  K.  Gerlach 

1992 

Phase  Codes 

Biphase  Perfect  Code 

S.W.  Golomb 

1992 

Biphase  Codes 

Ipatov 

V.P.  Ipatov 

1992 

Codes  with  minimal  peak  response  loss 

P(n,k) 

T.  Felhauer 

1994 

Phase  Codes 

Px  Code 

P.B.  Rapajic,  R.A.  Kennedy 

1998 

Phase  Codes 

PONS  based  Complementary  code 

P.  Zulch,  M.  Wicks,  et  al. 

2002 

Complementary  Codes 

Orthogonal  Codes 

N.  Levanon,  E.  Mozeson 

2003 

Train  of  Orthogonal  Coded  Pulses 

Multicarrier  Phase  Coded  Pulse 

N.  Levanon,  E.  Mozeson 

2002 

Multicarrier  Phase  Radar  Signals 

Table  VI 


Typical  codes  and  their  inventors. 

V.  Coherent  Train  of  Diverse  Pulses 

With  reference  to  a  train  of  coherent  pulses,  quite  often,  in  practice,  the  pulses  are  modu¬ 
lated  and  are  not  identical.  Modulation  produces  a  wider  bandwidth,  hence  pulse  compression. 
Moreover,  diversity  between  the  pulses  of  the  train  can  be  exploited  to  obtain  advantages  such 
as  lower  delay  sidelobes  or  lower  recurrent  lobes.  The  addition  of  a  modulation,  keeping  the 
pulses  identical,  leads  to  an  analytic  expression  for  the  ambiguity  function.  Adding  a  diversity 
in  amplitude  or  diverse  modulations  in  the  pulses  usually  requires  a  numerical  analysis  (but  for 
in  some  simple  cases),  since  the  ambiguity  function  is  only  available  as  an  integral  function. 
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Figure  10:  Partial  ambiguity  function  (|r|  <  tv)  of  a  coherent  uniform  train  of  eight  LFM  pulses  kti  =  20,  and 

Tr  =  9  tp. 


A  train  of  identical  LFM  pulses  provides  both  a  good  range  and  a  good  Doppler  resolution 
(probably  the  most  used  radar  signal  in  airborne  applications).  However  its  ambiguity  function 
(Figure  10)  still  presents  significant  sidelobes  both  in  delay  and  in  Doppler. 

Performing  an  interpulse  amplitude  weighting  permits  to  reduce  the  Doppler  sidelobes  at  the 
price  of  a  larger  Doppler  main  lobe  (at  v  =  0)  and  recurrent  lobes.  Similarly,  an  intrapulse 
weighting  in  LFM  mitigates  range  sidelobes.  The  aforementioned  weightings  can  be  thus  com¬ 
bined  to  reduce  both  range  and  Doppler  sidelobes. 

Another  method  for  diversifying  the  identical  pulse  train  relies  on  staggering  the  Pulse  Rep¬ 
etition  Frequency  (PRF)  obtaining  a  mitigation  of  the  blind  speeds  problem. 

Exploiting  a  pulse-to-pulse  diversity  can  lead  to: 

•  a  reduction  in  the  height  of  the  recurrent  (range)  lobes  of  the  autocorrelation  function  (i.e. 
around  r  =  nTR,  n  =  ±1,  ±2, . . .); 

•  a  reduction  of  near  range  sidelobes  (i.e.  around  |r|  <  tp); 

•  an  increase  in  the  overall  bandwidth  of  the  signal  while  maintaining  relatively  narrow 
instantaneous  bandwidth. 

In  this  context  we  mention  the  stepped-frequency  pulse  train  (which  consists  of  adding  a  fre¬ 
quency  step  A /  between  consecutive  pulses)  as  an  efficient  way  to  obtain  large  overall  bandwidth 
(hence  improved  range  resolution),  while  maintaining  relatively  narrow  instantaneous  bandwidth 
[36],  Moreover,  a  careful  selection  of  the  frequency  stepping  and  LFM  slopes  can  possibly 
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eliminate  the  ambiguous  peaks  of  the  autocorrelation  function  [37]-[38]  (of  course  complying 
with  the  unitary  volume  constraint). 


ORGANIZATION 


VI.  Conclusions 

In  this  lecture,  we  have  presented  an  overview  of  the  basic  theory  of  radar  waveforms.  First  of 
all,  we  have  provided  the  concept  of  ambiguity  function  and  have  discussed  its  relevant  properties. 
Then,  we  have  reviewed  the  most  common  radar  signals  together  with  their  ambiguity  functions. 
The  effects  of  a  possible  signal  coding  has  been  described:  a  special  attention  has  been  deserved 
to  the  class  of  frequency  coded  waveforms  through  a  Costas  sequence.  Their  construction,  relying 
upon  the  theory  of  Galois  fields,  as  well  as  the  desirable  properties  deriving  from  the  use  of 
a  Costas-pushing  sequence  are  discussed.  In  the  last  part  of  the  lecture,  we  have  considered 
coherent  pulse  trains  of  diverse  pulses  and  have  discussed  the  effects  of  modulating  the  different 
parameters  of  the  train. 

Before  concluding,  we  highlight  that  one  of  the  trade-offs  in  radar  signal  design  is  between 
constant  modulus  and  ambiguity  function  sidelobes.  Efficient  RF  power  amplifiers  are  presently 
operating  at  saturation  and  do  not  allow  linear  changes  in  amplitude.  On  the  other  hand,  sidelobe 
reduction  in  range  or  Doppler  usually  require  amplitude  variations  (weighting). 

Another  conflict  involving  a  linear  power  amplifier  relates  variable  amplitude  and  spectrum. 
The  spectrum  of  constant  amplitude  pulse  signal  exhibits  sidelobes  which  decay  very  slowly.  Such 
behavior  may  often  violate  spectrum  emission  regulations  and  can  possibly  cause  interference  to 
neighboring  radars  and/or  other  telecommunication  apparatus.  On  the  contrary,  suitable  variable 
amplitude  pulses  are  characterized  by  a  rapid  decay  of  the  spectral  tails. 

Further  radar  waveforms  with  variable  amplitude  are  the  Huffman-coded  signal  and  the  mul¬ 
ticarrier  signals  [1].  These  two  examples  show  that  removing  the  constant  amplitude  restriction 
provides  signals  with  additional  degrees  of  freedom  which  can  be  exploited  to  further  optimize 
the  system  performance. 
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